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1. ABSTRACT 


I This report addresses the development of a multidisciplinary optimization procedure using an 
efficient semi-analytical sensitivity analysis technique and multilevel decomposition for the design 
of aerospace vehicles. A semi-analytical sensitivity analysis procedure is developed for calculating 
computational grid sensitivities and aerodynamic design sensitivities. Accuracy and efficiency of 
the sensitivity analysis procedure is established through comparison of the results with those 
obtained using a finite difference technique. The developed sensitivity analysis techniques are then 
used within a multidisciplinary optimization procedure for designing aerospace vehicles. The 
optimization problem, with the integration of aerodynamics and structures, is decomposed into two 
levels. Optimization is performed for improved aerodynamic performance at the first level and 
improved structural performance at the second level. Aerodynamic analysis is performed by 
solving the three-dimensional parabolized Navier Stokes equations. A nonlinear programming 
technique and an approximate analysis procedure are used for optimization. The procedure 
developed is applied to design the wing of a high speed aircraft. Results obtained show significant 
improvements in the aircraft aerodynamic and structural performance when compared to a reference 
or baseline configuration. The use of the semi-analytical sensitivity technique provides significant 
computational savings. 

2. INTRODUCTION 

Analysis and design of aerospace vehicles are associated with complex multidisciplinary 
couplings. The development of an efficient optimization procedure for the design of aircraft must 
incorporate the interactions between disciplines such as aerodynamics, dynamics, aeroelastic 
stability, structures, controls and acoustics. However, the validity of the designs obtained using 
optimization techniques depends strongly upon the accuracy of the analysis procedures used and it 
is essential to integrate sufficiently comprehensive analysis procedures within the closed-loop 
optimization. Such procedures are computationally intensive, and therefore, can be prohibitive 
within an optimization environment. For example, it is essential to use a comprehensive 
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aerodynamic analysis procedure to solve the complex flow field associated with high speed 
aircraft. Over the past few years. Computational Fluid Dynamics (CFD) has evolved rapidly with 
the development of numerous numerical algorithms. Although accurate detailed analyses of many 
complex flow fields are now possible using supercomputers, viscous-compressible flow 
simulations of wing-body configurations can require several CPU hours per steady-state solution. 
Therefore, the use of such comprehensive analytical procedures for design optimization can be 
prohibitively expensive if a gradient-based technique is used. 

Sensitivity analysis, in which the derivative of a system performance function (e.g., the lift or 
drag of an aircraft) with respect to a design variable (e.g., a parameter controlling the wing 
planform) is calculated, is an essential ingredient in design optimization. A widely used technique 
for performing aerodynamic sensitivity analysis is the method of finite differences. The use of this 
method is associated with several calls to the flow analysis routine. Although this technique is 
conceptually simple, the associated computational cost is prohibitive when used in an optimization 
problem involving a large number of design variables. Therefore, it is necessary to develop 
efficient techniques to calculate aerodynamic sensitivities, so that advanced CFD codes may be 
more useful as practical design tools in multidisciplinary optimization environments. 

Two popular alternatives to the finite difference technique are the direct differentiation 
approach and the adjoint variable approach. These techniques are widely used in structural 
sensitivity calculations [1,2]. In both techniques, the actual governing equations are differentiated 
with respect to the design variables using chain rule. The direct differentiation approach yields a 
large system of equations involving the desired sensitivities that can be solved directly. In the 
adjoint variable approach, adjoint variables are obtained as the solution to an adjoint problem. The 
adjoint variables are then used to calculate the sensitivities. These two techniques are equivalent 
and yield identical results for the sensitivities. More recently, there has been widespread interest in 
using these techniques for calculating aerodynamic sensitivities. Carlson and Elbanna [3] have 
used the direct differentiation technique to differentiate the discretized transonic small perturbation 
equations and obtain aerodynamic sensitivities. Baysal et al. [4, 5] have performed discrete 
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sensitivity analysis using the Euler equations. Taylor et al. [6] and Newman et al. [7] have 
developed a semi-analytical sensitivity analysis procedure for the thin-layer Navier-Stokes 
equations using an incremental strategy. Jameson et al. [8] have proposed a continuous sensitivity 
approach using the adjoint variable method to calculate aerodynamic sensitivities. In the 
continuous sensitivity approach, the governing equations are differentiated prior to their 
discretization. The sensitivities are calculated using a numerical algorithm similar to the one used 
for obtaining the flow solution. Therefore, the continuous sensitivity approach needs to be 
modified, depending upon the governing equations that are differentiated. In the discrete 
sensitivity approach, the discretized governing algebraic equations are differentiated. Although 
there is a need for solving a large system of equations, this procedure can easily be adapted to 
different analysis procedures. In the present research, the aerodynamic sensitivities have been 
calculated by directly differentiating the discretized governing parabolized Navier Stokes (PNS) 
equations [9]. The discrete direct differentiation approach has been adopted over the continuous 
sensitivity approach because the CFD procedure used in this research to solve the PNS equations is 
based on a finite volume approach. This finite volume algorithm is more readily amenable to the 
direct differentiation approach than the continuous sensitivity approach. 

Two main ingredients in an aerodynamic sensitivity analysis procedure are: (1) the calculation 
of the sensitivities of the discretized flow variables and (2) the calculation of the sensitivities of the 
computational grid with respect to the aerodynamic design variables. It has been well recognized 
that the sensitivities of the flow variables are dependent upon the sensitivities of the computational 
grid [4-9]. However, in most of the aforementioned work, brute force finite difference techniques 
were used to calculate the grid sensitivities. Very few formal investigations have been reported on 
the development of analytical or semi-analytical techniques for computing grid sensitivities. High 
quality elliptic and hyperbolic grid generation codes are often used for generating meshes for 
aircraft configurations [10]. The use of the finite difference method for calculating grid 
sensitivities can be computationally prohibitive in such situations. Korivi et al. [11] developed a 
grid sensitivity analysis wherein the Jacobian matrix of the entire grid with respect to the grid 
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points on the boundary of the domain is calculated. The sensitivities of the surface grid points are 
calculated using an elastic membrane analogy to represent the computational domain, and the 
surface grid sensitivities are calculated from a structural analysis code using the finite elements 
method. Extension of this technique to complex three dimensional flow fields can be extremely 
complicated and time consuming. Further, the use of an additional structural analysis code 
increases computing time. Sadrehaghighi et al. [12] proposed an analytical approach for 
calculating grid sensitivities in which algebraic grid generation is performed using transfinite 
interpolation and surface parameterization in terms of design variables. The transfinite 
interpolation equations are analytically differentiated to obtain the grid sensitivities. The most 
general parameterization of the boundaries would require the specification of every grid point on 
the boundary. This, however, is impractical from a computational point of view. A quasi- 
analytical parameterization is used in Ref. 12 which allows the aircraft component to be specified 
by a relatively smaller number of parameters. However, the technique does not offer a great 
amount of generality because most CFD codes use complex grids which are generated using 
methods based on partial differential equations. 

In the present research, the grid sensitivity parameters are efficiently calculated, without any 
loss of generality and complexity, by directly differentiating the elliptic and hyperbolic grid 
generation equations [13]. This results in a large system of equations which can be solved readily 
to yield the grid sensitivities. The technique developed is not restricted to two dimensional 
problems and can be applied to three dimensional problems without any additional effort. The 
developed grid sensitivity technique is then used in conjunction with the semi-analytical 
aerodynamic sensitivity analysis procedure developed by Chattopadhyay and Pagaldipti [9], for 
calculating aerodynamic design sensitivities within a multidisciplinary optimization procedure. 

The development of the computationally efficient semi-analytical sensitivity analysis 
techniques allows the integration of comprehensive CFD-based codes within realistic 
multidisciplinary optimization studies. The necessity of multidisciplinary coupling in successful 
design optimization has been recognized. Recently, attempts have been made in coupling two or 
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more of these disciplines [14-18] in which optimization was performed by addressing all the 
design criteria in a single level. This “all-at-once” optimization procedure, in which all the 
disciplines are coupled inside a single loop and optimization is performed based on criteria 
involving every discipline, can be inefficient and time consuming. Decomposition techniques are 
often used to simplify such complex optimization problems into a number of sub-problems. 
Multilevel decomposition techniques have been applied to problems based on a single discipline 
[19-24] in structural applications. Recently, attempts have been made to use these techniques for 
multidisciplinary optimization of rotary wing aircraft. Adelman et al. [25] developed a two-level 
procedure for performing integrated aerodynamic, dynamic and structural optimization of rotor 
blades, based on the multilevel optimization strategy described in Ref. 22. Chattopadhyay et al. 
[26] developed a three-level procedure for optimization of helicopter rotor blades with the 
integration of aerodynamics, dynamics, aeroelastic stability, and structures. In a multidisciplinary 
design problem, the number of levels in a multilevel decomposition procedure typically depends 
upon the number of disciplines involved. Individual optimization is performed at each level using 
analysis procedures pertaining to that level. Optimal sensitivity parameters are exchanged between 
the levels to provide the necessary coupling. An optimal design is obtained when each individual 
level is converged and overall convergence is achieved. Therefore, the speed of obtaining a fully 
converged result depends upon the strength of coupling between the various levels. In the present 
research, the developed semi-analytical sensitivity techniques are used within a multilevel 
optimization procedure for designing aerospace vehicles with the coupling of aerodynamics and 
structures. 


3, OBJECTIVES 

The objectives of the current research project are as follows. 

(a) Development of a semi-analytical approach for calculating aerodynamic design sensitivities. 

(b) Development of a semi-analytical approach for calculating grid sensitivities to be used within 
the aerodynamic sensitivity analysis. 
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(c) Development of a multiobjective, multilevel optimization procedure for aerospace vehicles 
with the integration of the necessary disciplines such as aerodynamics and structures. 

These objectives have all been achieved successfully, details of which are furnished in the 
following sections. 

4. DISCRETE SEMI-ANALYTICAL AERODYNAMIC SENSITIVITY ANALYSIS 


The sensitivities of the aerodynamic coefficients of the aircraft with respect to its relevant 
geometric parameters are calculated in the present research using a discrete, direct differentiation 
approach. This approach is described in detail here. In general, an aerodynamic performance 
coefficient, Cj, depends on the steady-state flow variables, Q\ the vector of computational grid 
coordinates, X, and, sometimes, explicitly on the vector of independent design variables, <J>. 
Mathematically, 


Cj = Cj(Q*(<I>), X(d>), <D) 


(1) 


The derivative of Cj with respect to the i* design variable, 4> i5 is expressed as follows. 

*i= + tef® + acj 

d<t>i |dQ* | 1 3<)>i / |dX j 13<|>i j 3<|>i 


{ 3C 1 dC dC 

— U, — - and — — are easily calculated knowing the explicit dependence of C. 
dQ I dx d<pi J 

on Q*, X and (J),. The term [^-1, which represents the sensitivity of the steady state flow variables 

1 3<l>i I 


with respect to the i* design variable, is calculated using the direct differentiation technique. In the 
discrete sensitivity approach, the discretized flow equations are directly differentiated, as described 
next. The discretized flow equations which model the flow can be written as follows. 


(R(Q*(<D), X(<D),<D)} = {0} 


(3) 


Equation 3, differentiated with respect to <}>;, yields 



= p R i 

W! 


i T (axl 


d<j>ij 

|5Q‘l 

1 3<t>i | 

lax, 

( U<t>i / 
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( 4 ) 
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Equation 4 represents a set of linear algebraic equations in ^2- which can be solved easily. It is to 

4, can be calculated easily, knowing the explicit 


be noted that the terms 


— 1 and — in Eq. 

dQ ax I d<pi 


dependence of {R} on Q*, X and <J)j. 


5. GRID SENSITIVITY TECHNIQUE 

The term ||^j appearing in Eqs. 2 and 4 represents the grid sensitivity vector which must be 

computed semi-analytically. The semi-analytical grid sensitivity approach is illustrated on a 
hyperbolic grid generator. In general, a three dimensional hyperbolic grid generation code 
generates a two-dimensional grid at various stations along the longitudinal direction by solving the 
following equations [10]. 


dydy [ dzdz _p 

ana c + 

dydz dzdy _ r 


(5) 


( 6 ) 


where the xyz coordinate system is a Cartesian coordinate system fitted to the body and £nC 

coordinate system is the computational domain used by the CFD procedure for aerodynamic 

analysis [Fig. 1]. In Eq. 6, F(n,C) is a known function approximating the Jacobian of 

transformation between the xyz and the £nC coordinate systems. Equations 5 and 6 are discretized 
and solved numerically to obtain the grid vector X. The grid sensitivity vector, (— ], can be 


d<t>i 

obtained by directly differentiating Eqs. 5 and 6 with respect to <|> t after their discretization, as 
follows. 


fe-U d yj,k \ / dyj,k+l dyjk\ 

dy \ dfo d(f)j ) | 3y\ d<[>i d<}>i / _ dz 


f dz j+l,k dz j,k ^ 


f dz j,k+l dzjkl 


i 


#i / + 3z\ d( )>i d<>i 


^ fe+i-Cj) ac (^ + i-qk) an (; j+1 -g) 3^ (%+i-nk) 

/ dz j+l,k dz jkl 


= 0 (7) 


I d yj,k+i _ d yj,k \ I d yj+u d yj,k \ l dz j,k+i dz j,k 

3y \ d(j>j d<j)j j | 8z\ dfo d<j)j ) 3z\ d<j>j d(>i / 3y\ dfo dfo _ dF 

(Cj+i - Cj) a? (m+i-nk) 'an (Cj+i - Cj) ac 
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Equations 7 and 8 represent a system of equations which can be solved readily to yield the grid 
sensitivity vector, 

The grid sensitivity and the aerodynamic sensitivity techniques yield large systems of algebraic 
equations characterized by sparse coefficient matrices. Direct methods for solving these systems of 
equations are inefficient and iterative techniques have been used in this work. The systems of 
equations are solved using the successive over relaxation (SOR) scheme. A detailed description of 
this iterative technique can be found in Ref. 10 and is not repeated here. 

6. PROBLEM FORMULATION USING MULTILEVEL DECOMPOSITION 

This section describes the multilevel decomposition technique useful in breaking up coupled 
and complex optimization problems into simpler optimization problems. The multilevel 
decomposition procedure is illustrated through a two-level formulation. Each level is a 
multiobjective optimization problem characterized by a vector of objective functions, constraints 
and design variables. During optimization at a particular level, it is essential to maintain the 
objective functions and design variables of lower levels close to their optimum values. Therefore, 
constraints are imposed on the perturbations to the lower level objective functions and design 
variables to prevent significant changes. These parameters are called optimal sensitivity 
derivatives, and they establish the necessary link between the various levels of optimization. The 
multilevel decomposition procedure is outlined below. 

Level 1: 

Minimize 

FjOJ) 1 ) i = 1, ..., NOBJ 1 

subject to 

gk^ 1 ) < 0 k = 1, ..., NC 1 
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j = 1, NOBJ 2 


•TOV 1 

X 






* e 2j 


1L 1 1U 

<t>i ^ <l>i ^ <t>i 


i = 1, NDV 1 




NDV 2 


where F l and F 2 are the objective function vectors at levels 1 and 2 respectively, g 1 and g 2 are the 
corresponding constraint vectors and (J) 1 and (j ) 2 are the corresponding design variable vectors. The 
quantity £2j is a tolerance on the change in the j 111 objective of level 2 during optimization at level 1. 
Superscripts L and U represent lower and upper bounds respectively, and superscript * represents 


optimum values obtained at level 2. 


Finally, the quantities. 




and 


a£ 

d<J>i 


are the optimal 


sensitivity parameters of level 2 objective function and design variable vectors, respectively, with 
respect to the level 1 design variables. 


Level 2: 


Minimize 

FjCtt) 1 *, <|> 2 ) i = 1..NOBJ 2 

subject to 

gjfo) 1 *, <l> 2 ) ^ 0 k = 1..NC 2 

<t> 2L < <(>? < <|)f U i = 1..NDV 2 

where (j) 1 * is the optimum design variable vector from level 1. This vector is kept fixed during 
optimization at level 2. The optimization procedure cycles through the two levels before global 
convergence is achieved. A "cycle" is defined as one complete sweep through the two levels of 
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optimization. Optimization at an individual level also requires several "iterations" before local 
convergence is achieved. Cycling between the two levels is necessary to account for the coupling 
between the objective functions, constraints, and design variables pertaining to the levels. 

7. ANALYSIS 

The parabolized Navier-Stokes equations (PNS equations) have been used for the evaluation 
of three-dimensional, supersonic, viscous flow fields. The PNS equations are obtained from the 
full Navier-Stokes equations based on the following assumptions: (a) steady state, (b) the 
streamwise viscous gradients are neglected, and (c) the streamwise pressure gradient in the 
subsonic portion of the viscous flow near the body surface are approximated. The inviscid region 
of the flow field must be supersonic and the streamwise velocity component must be positive 
everywhere. Thus streamwise flow separation is not allowed but crossflow separation is allowed. 
Efficiency in computational time and memory requirements are achieved because the equations can 
be solved using a space-marching technique. The computational procedure used in this study, as 
implemented in the code, UPS3D [27] integrates the PNS equations using an implicit, 
approximately factored, finite-volume algorithm where the crossflow inviscid fluxes are evaluated 
by Roe’s flux-difference splitting scheme [28]. The UPS 3D code also has the capability to 
calculate the inviscid flow field by solving the PNS equations without the viscous terms. The 
upwind algorithm is used to improve the resolution of the shock waves over that obtained with the 
conventional central differencing schemes. A hyperbolic computational grid is used with 75 grid 
points along the circumferential (rj) direction, 80 grid points along the normal (Q direction, and 31 
grid points along the longitudinal (£) direction. Further refinement of the grid does not change the 
flow solution significantly. 

The aircraft wing structural analysis is performed using an inhouse code. The code is capable 
of analyzing multicelled box beams of arbitrary cross-section and tapered planform. The wing 
section is represented by a diamond airfoil (Fig. 2). The principal load carrying member in the 
wing is modeled by an isotropic box beam with a rectangular cross-section and unequal wall 
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thicknesses. The wall thicknesses are represented as fractions of local chord (c) as shown in 
Figure 2. The aircraft weight (W) is calculated as the sum of the weight of the box beam (Wbox) 
and the weight of the skin (W s ian)- The stresses (o) are calculated using thin wall theory. 

8. OPTIMIZATION 


A gradient based optimization technique based on the method of feasible directions [29] has 
been used to solve the optimization problems at levels 1 and 2. Since the optimization process 
requires several evaluations of the objective function and the constraints before an optimum design 
is obtained, the process can be very expensive if actual analyses are performed for each function 
evaluation. The objective function and constraints at levels 1 and 2 are, therefore, approximated 
using a two-point exponential approximation [30] based on the design variable values from the 
optimizer and the sensitivities of these functions for the current and the previous design points. 
This approximation is given by the following expression. 

a ND\ 

= mi) + i 

n=l 



- 1.0 


<l>ln ^F 


Pn 


«E»i) 


(9) 


where ^(O) is the approximation of the function F(d>) in the neighborhood of the current design 
variable vector, <3>i. The quantity (J) n is the n* design variable from the design variable vector <t>. 
NDV is the total number of design variables. The approximate values for the constraints, gj(<t>), 
are similarly calculated. The exponent p n is defined as 


Pn = 



+ 1.0 


( 10 ) 


where d>j refers to the design variable vector from the current cycle and Oq denotes the design 
variable vector from the previous cycle. Equation 9 indicates that in the limiting case of p n = 1, the 
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expansion is identical to the traditional first order Taylor series and when p n = -1, the two-point 
exponential approximation reduces to the reciprocal expansion form. Therefore, the exponent (p n ) 
can be interpreted as a “goodness of fit” parameter which explicitly determines the trade-offs 
between traditional and reciprocal Taylor series based expansions. p n is chosen to be within the 
interval, - 1 < p n < 1 thus resulting in a hybrid approximation technique. 

9. APPLICATION 

A two level optimization procedure has been developed for simultaneous improvement of 
aerodynamic and structural characteristics of high speed aircraft. From aerodynamics point of 
view, the aircraft drag coefficient (Cd) must be minimized while maintaining the lift coefficient 
(Cl) above a desired level. From a structural perspective, it is desirable to minimize the aircraft 
weight (W) while maintaining the stresses in the load carrying members within material limits. A 
semi-analytical technique is developed and used within the optimization procedure for calculating 
sensitivities of aerodynamic coefficients and the computational grid used in the aerodynamic 
analysis. The developed two level optimization procedure and sensitivity analysis technique are 
applied to a delta wing-body configuration related to high speed aircraft. 

9.1 Aircraft Configuration 

The two level optimization procedure and the semi- analytical sensitivity analysis techniques 
are applied to the design optimization of a delta wing-body configuration illustrated in Figure 2. 
This center body is characterized by a nose region where the radius increases parabolically from 
zero to its maximum value (r m ). Beyond the point where maximum radius is reached, the body 
remains a cylinder through the length of the aircraft. The delta wing is characterized by a leading 
edge sweep (X), root chord (c 0 ) and a wing span (w s ). The wing cross-section is a diamond airfoil 
[Fig. 1] with a thickness-to-chord ratio denoted tc. The wing leading edge sweep, root chord, 
wing span and thickness-to-chord ratio are used as design variables within the aerodynamic 
optimization at level 1. The load carrying structural member of the wing is modeled as a single 
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celled, isotropic box beam [Fig. 3]. The beam width-to-chord ratio (w c ), the horizontal wall 
thickness-to-chord ratio (ti) and the vertical wall thickness-to-chord ratio (t2) are used as design 
variables during the structural optimization at level 2. 

9.2 Two Level Optimization Problem 

In a typical aircraft design, the aerodynamic performance criteria dictate the planform shape of 
the aircraft wing. After designing the planform, the wing load carrying member is designed to 
carry the loads arising from aerodynamics, inertia and gravity. Keeping this in mind, the two level 
optimization procedure improves the aircraft aerodynamic performance at level 1 and the structural 
performance at level 2. These optimization problems are stated as follows. 

Level 1 ( Aerodynamic sV. 

Minimize 

Drag coefficient (Cd) 
subject to the constraints 

Cl ^ CL m j n (Lift coefficient constraint) 

W < W* (Optimal weight constraint) 

* 

L * NDV1 3 ( t )- 9 U 

<t> j2 * <^2 + s r-^ii * <t> j2 

1 ®Tii j = L.NDV2 (Structural design variables constraint) 

L u 

^ii - ^ii - ^ii i = l..NDVi (Side constraints on aerodynamic design variables) 

where <3>j denotes the vector of aerodynamic design variables, <J> 2 denotes the structural design 
variable vector, NDVi and NDV2 denote the number of aerodynamic and structural design 
variables respectively and the superscripts L and U denote lower and upper bounds on the design 
variables respectively. 
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Level 2 (Structures); 

Minimize 

Weight (W) 

subject to the constraints 

o < o all (Stress constraints on wing load carrying member) 

L U 

- ^i 2 - ^i2 i=l..NDV 2 (Side constraints on structural design variables) 


where o is the vector of stresses at specific locations on the aircraft wing and a a n is the 
corresponding vector of material limits. During the structural optimization at the second level, the 
aerodynamic design variables of level 1 (<()i) are held fixed at their optimum values. The quantities, 


aw* 

a^ii 


and 



are the sensitivities of the optimum values of the second level (structural) objective 


function and the design variables with respect to first level (aerodynamic) design variables. These 
optimal sensitivity derivatives provide the necessary coupling between the two disciplines. 


9.3 Results and Discussion 


The two level optimization procedure and the semi-analytical sensitivity analysis techniques 
are applied to the design optimization of the delta wing-body configuration [Fig. 2]. Results from 
the grid sensitivity technique and the aerodynamic sensitivity analysis are presented first followed 
by results obtained from the two level optimization. The grid sensitivity technique is applied to the 
three-dimensional hyperbolic grid around the delta wing-body configuration. The values of the 
design variables used are: wing root chord (c 0 ) = 7.08 m, leading edge sweep (A.) = 66.0 degrees, 
wing span (w s ) = 2.96 m and the wing thickness-to-chord ratio (t c ) = 0.052. The computational 
grid includes 75 grid points in the circumferential (n) direction, 31 grid points in the longitudinal 
(£) direction and 80 grid points in the normal (Q direction. The grid sensitivities obtained from the 
direct differentiation technique and finite difference technique are presented in Table 1. As shown, 
there is excellent agreement between the two techniques. Further, a comparison of the CPU time 
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shows a 42 percent reduction achieved for one complete grid sensitivity analysis using the 
developed procedure [Fig. 4]. This clearly demonstrates the significant computational savings 
achievable by using this approach, especially in a formal design optimization procedure where 
several such design sensitivity analysis are necessary. 

Results obtained from the aerodynamic sensitivity analysis procedure are presented next The 
analysis is performed for an angle of attack of 5 degrees and a flight Mach number of 2.5. The 
sensitivities of the drag coefficient (C D ) and the lift coefficient (Cl), calculated using the direct 
differentiation technique as well as the finite difference technique, are presented in Tables 2 and 3 
respectively. It must be noted that column 3 in Tables 2 and 3 presents the results of the semi- 
analytical aerodynamic sensitivity approach with finite difference grid sensitivity while column 4 
presents the results of the semi-analytical aerodynamic sensitivity approach with semi-analytical 
grid sensitivity. As shown, the results from both techniques are in excellent agreement. For one 
complete sensitivity analysis, the direct differentiation technique with finite difference grid 
sensitivity calculations results in a 30 percent reduction in computing time from the finite difference 
technique [Fig. 5]. The semi-analytical sensitivity analysis technique with semi-analytical grid 
sensitivity calculations yields a 45 percent reduction in computing time from the finite difference 
approach [Fig. 5]. This further illustrates the efficiency of the discrete semi-analytical technique 
for grid sensitivity calculations. In order to investigate the trend in CPU savings, the results are 
plotted for varying number of design variables. Figure 6 compares the computing time required, 
for one complete sensitivity analysis, using the finite difference technique, the semi-analytical 
aerodynamic sensitivity analysis using finite difference grid sensitivity, and the semi-analytical 
aerodynamic sensitivity analysis using semi-analytical grid sensitivity, as a function of the number 
of design variables. It is observed that, although all three curves follow a near linear relationship 
with the number of design variables, their slopes are quite different. The design sensitivity 
analysis with semi-analytical grid sensitivity calculations results in the minimum slope. This 
indicates that, although, as expected, the computation time increases almost linearly with the 
increase in the number of design variables, the increase is less significant in the semi-analytical 
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approaches and least in the case where analytical grid sensitivities are used. This further illustrates 
the efficiency of the grid sensitivity analysis technique. 

The results obtained from the two level optimization procedure are presented next. The 
reference values for the aerodynamic design variables are the same as those used to demonstrate the 
grid and aerodynamic sensitivity analysis. The load carrying structural member is modeled as an 
isotropic box beam made of 2014-T6 Aluminum alloy. The reference values of structural design 
variables are: spar width-to-chord ratio (w c ) = 0.5, spar horizontal wall thickness-to-chord ratio 
(ti) = 0.0015 and vertical wall thickness-to-chord ratio (t 2 ) = 0.0075. The iteration histories of the 
drag coefficient and the lift coefficient of the aircraft, during the level 1 aerodynamic optimization, 
are presented in Figs. 7 and 8 respectively. Significant improvements are observed in both the 
quantities. The drag coefficient decreases by 5.5 percent and the lift coefficient increases by 5.43 
percent. The lift coefficient is maintained at its improved value throughout the optimization 
procedure because, the lift coefficient constraint, imposed during the aerodynamic optimization, 
remains an active constraint. The optimizer, based on the method of feasible directions, remains 
on this constraint boundary, as it improves the objective function value. The constraint imposed 
on the aircraft weight, at level 1, is well satisfied during the level 1 optimization. Table 4 compares 
the reference and the optimum values of the aerodynamic design variables used in the aerodynamic 
optimization. The root chord is increased significantly from its reference value (8.76 percent), 
whereas the wing thickness-to-chord ratio is decreased significantly (40.29 percent). The wing 
span and the leading edge sweep are maintained close to their reference values (2.94 percent and 
1.01 percent respectively). The reduction in the wing thickness-to-chord ratio decreases the form 
drag of the diamond airfoil section. The optimization is driven by this reduction in drag. The 
increase in the wing planform area, caused by the increase in the root chord and the wing span, 
helps improve the lift of the aircraft in spite of the decrease in the wing thickness-to-chord ratio. 

The iteration history of the aircraft weight is presented in Fig. 9. Considerable reduction 
(18.13 percent) in the aircraft weight is observed from the reference to the optimum. The shear 
and the normal stresses at the blade root remain well within the allowable limits of the chosen 


16 



Aluminum alloy. Table 5 compares the reference and optimum values of the structural design 
variables used in the structural optimization. The optimum spar width to chord ratio, horizontal 
and vertical wall thickness to chord ratios decrease significantly from their reference values (60 
percent, 60 percent and 65.2 percent respectively). As mentioned above, the optimum 
configuration has an increased planform area compared to the reference configuration. It must be 
noted that the optimum weight is still lower than the reference due to the significant reduction in the 
wall thicknesses. 


10. SIGNIFICANCE OF THE RESEARCH 

An efficient, semi-analytical sensitivity analysis technique, based on the direct differentiation 
approach, has been developed for calculating aerodynamic and grid sensitivities. A 
multidisciplinary optimization procedure for designing wing-body configurations has been 
developed based on a two level decomposition. The aircraft aerodynamic and structural design 
requirements have been coupled within this two level optimization procedure. The following 
observations are made. 

1 . The results from the semi-analytical grid sensitivity and aerodynamic sensitivity analysis 
techniques compare very well with those obtained using a finite difference approach. 

2. The developed sensitivity analysis techniques yield significant savings in computational time 
and allow the use of comprehensive CFD within gradient based optimization procedures. 

3 . The semi-analytical aerodynamic sensitivity technique in conjunction with the developed grid 
sensitivity procedure shows the minimum increase in CPU time with increase in total number 
of design variables. 

4. The two level decomposition optimization procedure yields significant reductions in the 
aircraft drag and weight while improving the lift. 

5 . The reduction in the aircraft drag is predominantly due to the significant reduction in the wing 
thickness-to-chord ratio. The improvement in the lift is due to the increased planform area 
caused by the increases in the wing root chord and wing span. 
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Table 1. Grid sensitivity of the three-dimensional hyperbolic grid 


Grid point (x, y, z) 

Design variable 

Finite difference grid 
sensitivity method 

Direct differentiation 
grid sensitivity method 

(0.300, 0.044, 0.008) 

(16.406, 13.297, 
2.753) 

Sweep (A) 

Root chord (c 0 ) 
Wing span (w s ) 
Thickness/chord (t^ 

Sweep (A) 

Root chord (c 0 ) 
Wing span (w s ) 
Thickness/chord (t*.) 

(0.0, 0.0, 0.0) 

(0.0, 0.0, 0.0) 
(0.0, 0.0, 0.0) 
(0.0, 0.0, 0.0) 
(0.0, 0.0338, 0.3930) 

(0.0, 0.2886, 0.3356) 
(0.0, 0.8014, 0.9320) 
(0.0, 46.640, 54.240) 

(0.0, 0.0, 0.0) 

(0.0, 0.0, 0.0) 
(0.0, 0.0, 0.0) 
(0.0, 0.0, 0.0) 
(0.0, 0.0338, 0.3930) 

(0.0, 0.2884, 0.3355) 
(0.0, 0.8013, 0.9319) 
(0.0, 46.550, 54.132) 

Table 2. Sensitivity of the drag coefficient, (dCjydfy) 

Design variable 

Finite Difference 

Semi-analytical method 
(with finite difference 
grid sensitivity) 

Semi-analytical method 
(with semi-analytical 
grid sensitivity) 

Sweep (A) 
Root chord (c Q ) 
Wing span (w s ) 
Thickness/chord (t c ) 

-0.0097712 

0.0356022 

0.0239097 

1.4250140 

-0.0106882 

0.0353581 

0.0210748 

1.4219202 

-0.0106453 

0.0353786 

0.0229120 

1.4820189 


Table 3. Sensitivity of the lift coefficient, (dO/dt));) 


Design variable 

Finite difference 

Semi-analytical method 
(with finite difference 
grid sensitivity) 

Semi-analytical method 
(with semi-analytical 
grid sensitivity) 

Sweep (A) 

-0.0943048 

-0.0920100 

-0.0880802 

Root chord (c 0 ) 

0.0775282 

0.0727289 

0.0755285 

Wing span (w s ) 

0.0695290 

0.0628465 

0.0647937 

Thickness/chord (t c ) 

-1.8981815 

-1.7628147 

-1.8495492 
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Table 4. Aerodynamic design variables 


Design variable 

Reference 

Optimum 

Sweep (A.) 

66.00 deg. 

65.33 deg. 

Root Chord (c 0 ) 

7.08 m 

7.70 m 

Wing Span (wg) 

2.96 m 

3.047 m 

Thickness-to-chord ratio (tc) 

0.05200 

0.03105 


Table 5. Structural design variables 

Design variable Reference 

Optimum 

Spar width-to-chord ratio (w c ) 0. 5 

0.2 

Horizontal wall thickness-to-chord 0.0015 

ratio (tj) 

0.0006 

Vertical wall thickness-to-chord 0.0075 

ratio (t 2 ) 

0.00261 
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Nose length (In) 



c 


Figure 1. Coordinate systems 



Diamond airfoil (Section AA) 



Wing starting location (xw) 


Figure 2. Delta wing-body configuration 
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Wing span (Ws) 




Wing Spar (Box Beam) 


Figure 3. Wing cross-section and wing spar (box beam) 


B Direct differentiation 
■ Finite difference 



Figure 4. Comparison of CPU time for the sensitivity of hyperbolic grid 
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CPU time required CPU time on CRAY-YMP 


□ Finite difference 

E3 Semi- Analytical Method with Finite Difference Grid Sensitivity 
■ Semi- Analytical Method with Semi- Analytical Grid Sensitivity 



Figure 5. Comparison of CPU time for aerodynamic sensitivity analysis 


Finite difference 



Figure 6. Comparison of CPU time with design variables 
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Iteration 


Figure 8. Iteration history of lift coefficient 
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Figure 9. Iteration history of aircraft weight 
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